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ABSTRACT 

Through  the  use  of  asymptotic  and  perturbation  methods, 
this  thesis  presents  a  theoretical  study  of  the  flow  of  a 
stratified  fluid  over  variable  topography  as  a  model  of  the 
resonant  interaction  of  a  submarine's  wake  with  a  stratified 
fluid.  Such  resonant  interactions  may  be  able  to  produce 
significant  upstream  disturbances.  The  long  time  solution 
obtained  in  our  model  exhibits  growth  in  time  for  the  resonant 
case,  indicating  that  perhaps  nonlinear  effects,  balanced  by 
dispersion,  could  cause  significant  upstream  disturbances. 
These  disturbances  could  allow  the  possibility  of  detecting  a 
submarine  by  the  internal  waves  its  wake  generates  in  the 
resonant  case. 
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I .   INTRODUCTION 

The  nonacoustic  detection  of  submarines  is  an  area  of 
intensive  and  well  funded  research.  Although,  it  is  not  clear 
whether  devices  based  on  nonacoustic  means  of  detection  will 
become  as  useful  as  those  based  on  acoustic  means,  the 
military  importance  of  such  devices  demands  pursuing  such 
avenues  of  technological  research.  One  possible  nonacoustic 
means  of  detection  is  to  attempt  to  detect  the  submarine 
optically  through  the  use  of  a  laser  by  examining  the  way  the 
submarine  reflects  or  absorbs  blue-green  light.  Another  more 
promising  method  is  to  try  to  detect  the  submarine  by  the  heat 
it  emits.  However,  an  even  more  promising  technique  is  the 
detection  of  submarines  by  the  internal  waves  they  generate. 

In  the  ocean  a  stable  density  stratification  generally 
exists.  As  the  submarine  moves  through  the  density-stratified 
ocean,  it  leaves  behind  a  wake  capable  of  generating  internal 
waves.  These  internal  waves  result  from  the  stabilizing 
effect  of  the  stratification:  as  a  water  particle  is 
displaced,  it  experiences  a  buoyant  restoring  force  in  the 
opposite  direction  from  which  it  was  displaced. 

Many  theoretical  and  experimental  studies  have  been 
conducted  on  the  subject  of  the  flow  of  a  stratified  fluid 
over  a  localized  topography.  For  the  two-dimensional  flow  of 
an  inviscid,  incompressible  stratified  fluid,  most  studies 


have  been  confined  to  steady  flow  or  linearized  theory. 
Earlier  studies  concerned  themselves  only  with  a  description 
of  the  downstream  stationary  lee-wave  field.  A  question  that 
has  generated  some  controversy  is  the  phenomenon  of  upstream 
influence.  To  what  extent  do  disturbances  generated  by  the 
topography  through  transient  and  nonlinear  processes  propagate 
upstream? 

Mclntyre  (19  72)  showed  that  the  only  significant  upstream 
disturbances  were  weak  second  order  long -wave  motions 
generated  by  nonlinear  interactions  in  the  lee- wave.  Mclntyre 
(1972)  extended  the  earlier  theories  of  Benjamin  (1970)  and 
Keady  (19  71)  by  performing  an  expansion  using  as  the  small 
parameter  the  ratio  of  the  height  of  the  topography  to  the 
total  fluid  depth.  Baines  (1977,1979),  on  the  other  hand,  in 
experimentally  examining  the  flow  of  a  continuously  stratified 
fluid  over  topography,  found  upstream  disturbances  that  were 
first  order  in  the  topography  height.  These  upstream 
disturbances  were  not  only  generated  by  nonlinear  processes 
over  the  topography  but  also  travelled  at  the  long -wave  speed. 
Baines  also  observed  that  these  upstream  disturbances  were 
very  strong  when  the  fluid  approached  resonance.  Here, 
resonance  is  defined  as  the  near  coincidence  of  the  basic  flow 
speed  and  one  of  the  free  long- wave  speeds.  Grimshaw  and 
Smyth  (1986)  presented  a  theoretical  study  of  stratified  fluid 
flow  over  localized  topography  for  the  two-dimensional  case 
when  the  basic  flow  is  near  resonance.    One  interesting 


feature  of  their  results  is  the  appearance  of  significant 
upstream  disturbances  as  observed  by  Baines  (1977,1979). 
Grimshaw  and  Smyth  showed  that  solitary  waves  were  formed  when 
the  localized  forcing  over  the  topography  produced  a  positive 
disturbance  upstream,  and  an  oscillatory  wavetrain  was  formed 
when  the  disturbance  had  negative  polarity. 

Motivated  by  these  results,  this  thesis  presents  a 
theoretical  study  of  the  flow  of  a  stratified  fluid  over 
localized  topography  for  the  three  dimensional  case. 


II.  PROBLEM  FORMULATION 

A.   DIMENSIONAL  PROBLEM 

We  consider  the  three-dimensional  flow  of  an  inviscid, 
incompressible  stratified  fluid.  Figure  1  (see  page  5) 
provides  a  sketch  of  the  coordinate  system.  We  assume  that 
the  basic  state  has  a  constant  horizontal  velocity,  V,  from 
left  to  right  with  equilibrium  pressure  and  density 
respectively  as  p0(z)  and  p0(z).  Here  and  throughout  this 
paper  barred  variables  represent  dimensional  variables.  We 
introduce  f  as  the  vertical  particle  displacement  from  some 
initial  position  at  time  t  =  0  to  the  position  at  time  t  =  T. 
Using  the  Lagrangian  representation  f(x,y,z,T)  =  z^  -  z0,  the 
density,  p(x,y,z,t),  becomes  p0(z- f(x,y, z, t) ) .  Therefore,  the 
equations  of  motion  are 


u-  +  v-  +  w-z  =   0,  (2.1) 


po(z-n^    +  P-X   =  0,  (2.2 


P<,(z-f)£r  +  p-y  =  0,  (2.3) 


p„u-n|^  +  P„u-r>sr  +  pe  =  o,  (2.4 


1 

z 

=  H  +  t/ 

V 

1 
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/    y 

/z  =  f                      / 
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H 

► 

Figure  1.  The  Coordinate  System 


w  =    El,  (2.5) 


where 


_£  =  A.    +  (V+u)-4  +  V-4  +  w— .  (2.6) 

DT    dT        dx     dy    dz 


Here,  u,  v,  and  w  are  the  velocity  components  relative  to  the 
basic  state  which  has  velocity  (V,0,0)  with  respect  to  a 
right-handed  Cartesian  coordinate  system. 

The  bottom  topography  is  given  by  z  =  f(x,y,t),  where 
initially  f(x;y,0)  =  0.  For  large  times  the  topography 
attains  a  steady  value  g(x,y)  and  hence  we  assume  that 

f(x,y,t)   ~g(x,y)   ast-*oo.  (2.7) 

We  also  assume  that  at  all  times  the  topography  is  localized 
where  f  (x,y,t)  ^  0  as  x2+y2  -*  oo.  The  kinematic  bottom  boundary 
condition  is  then 

7  =  ?    at    ~z  =1.  (2.8) 


We  are  assuming  that  this  "topography"  is  created  by  the 
submarine's  wake,  which  when  viewed  from  the  moving  submarine, 
actually  has  a  steady  shape.  In  this  thesis,  we  investigate 
the  interaction  of  the  wake  with  internal  waves.  Kallen 
(1987)  has  investigated  this  type  of  interaction  in  a  linear, 
steady  state  model  and  comments  in  detail  on  the  validity  of 
assuming  that  the  wake  acts  as  rigid  topography.  Furthermore, 
from  experimental  results,  Lin  and  Pao  (1979)  show  the  mean 
flow  occurs  at  the  wake  boundaries  implying  that  part  of  the 
flow  travels  around  rather  than  through  the  wake.  It  is  this 
flow  around  the  wake  that  generates  the  internal  waves.  Our 
model,  in  fact,  does  not  depend  critically  on  the  assumed  form 
of  "topography"  the  submarine  creates  and  our  purpose  here  is 
to  investigate  a  time  dependent  model. 

The  upper  boundary,  assumed  to  be  a  free  surface,  is 
located  at  z  =  H  +  r?(x,y,t),  where  H  is  the  vertical  distance 
from  the  surface  to  the  submarine  and  ry(x,y,t)  is  the  free 
surface  displacement.  One  must  remember  to  include  not  only 
the  kinematic  boundary  condition  at  the  upper  boundary  but  the 
dynamic  boundary  condition  as  well.  The  upper  boundary 
conditions  take  the  form: 

£  =  V  at    z  =  H  +  ~rj,  (2.9a) 

po(H+v)    +  PJH+v)    =P0(H)       at      ~^  =  H  +  ^.  (2.9b) 


Here,  pp  is  the  perturbation  to  the  pressure  p0  of  the  basic 
state. 

B.   DIMENSIONAL  SCALING 

The  solution  of  practical  problems  in  fluid  mechanics 
requires  both  theoretical  developments  and  experimental 
results.  By  grouping  significant  dimensional  quantities  into 
dimensionless  parameters  it  is  possible  to  present  results  in 
a  compact  form  which  is  applicable  to  all  similar  situations. 

In  the  problem,  the  characteristic  horizontal  length  scale 
L  is  determined  by  the  ratio  U/N,  where  U  is  a  characteristic 
horizontal  velocity  and  N  is  the  Brunt- Vaisala  frequency  at 
the  location  of  the  submarine.  We  shall  take  the 
characteristic  horizontal  velocity  U  to  be  given  by  N,H,  (U  = 
N,H)  ,  where  N,  is  the  maximum  Brunt -Vaisala  frequency  and  H, 
the  characteristic  vertical  length,  is  the  depth  of  the 
submarine.  Typically  N,  is  10"2  sec"1  and  H  is  200  meters  (m)  , 
so  U  is  2  m/sec.  This  is  a  typical  order  of  magnitude 
patrolling  speed  for  a  submarine  at  that  depth.  Of  course, 
0(1)  variations  in  N,  and  H  can  be  made  without  changing  our 
analysis . 

From  the  mass  conservation  equation  and  with  the  imposed 
horizontal  length  and  velocity  scales,  we  see  that  the 
vertical  velocity  scale,  W,  is  UH/L.  The  horizontal  length 
and  velocity  scales  are  also  used  in  determining  the 
characteristic  time  scale,  T,  where  T  is  defined  by  T  =  L/U  = 


1/N.  Lastly,  the  characteristic  pressure  scale  is  defined  as 
p  =  /o,gH  where  p,  is  a  characteristic  density.  These  scales 
define  the  parameter  (3  =  N,2H/g,  which  is  a  dimensionless 
measure  of  the  ratio  of  the  buoyancy  force  to  the 
gravitational  force. 

C .   NOND IMENS I ONAL I Z AT I ON 

By  using  the  various  dimensional  scales  introduced  in 
section  B,  we  shall  nondimensionalize  the  equations  of  motion 
(2.1-2.6)  along  with  the  boundary  conditions  (2. 8-2. 9b). 
However,  before  deriving  the  nondimensional  equations  of 
motion  and  boundary  conditions,  we  will  consider  several 
related  dimensional  equations  that  will  assist  in  simplifying 
the  nondimensional  derivation  of  equations  (2.1-2.6)  and 
boundary  conditions  (2. 8-2. 9b).  We  first  introduce  the 
following  nondimensional  variables: 


p  -  p° 

H 

Po    -- 

Po 

Po      P,gtf 

P\ 

n-  I    , 

f  -1  . 

3          H 

(2.10) 


To  nondimensionalize  the  hydrostatic  equation,  we  begin  with 


dp 

-£    =  ~9Po-  (2.1D 

dz 


Using  (2.10)  yields 


Po(z)    =  -Po.  (2.12) 


The  dimensional  Brunt -Vaisala  frequency  is  defined  as 


g  df>o 


N'  =   -JL-LS  (2.13) 

Po  dz 


where  dp0/dz  <  0.   Using  (2.10),  we  obtain 


N2H   ,     x  dp0 

-J—N2  =   -—-LI  (2.14 

9  Po  dz 


but  /3  =  N,2H/g,  therefore  (2.14)  simplifies  to 


T2   _       dPo 


Pof3N2  =  --11.  (2.15) 

dz 
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In  nondimensionalizing  the  Lagrangian  representation  of  the 
particle's   vertical    displacement,    we   start   with 

~zl}  =  z  -  1  (x,y,z,t)  .  (2.16) 

Using    (2.10)    yields 

z0  =  z  -   {(x,y,z,  t)  .  (2.17) 

The  last  equation  needed  before  considering  the  equations  of 
motion  and  boundary  conditions  is  the  dimensional  pressure 
equation: 

p(x,y,~z,t)    =  pxgHp0{z)    +   p1gHj3g(x;y/  z,  t)  (2.18) 

where  pp  =  p,gH/Jq  (x,y,  z,  t)  and  0  is  defined  after  (2.14). 
Thus,  q(x,y,z,t)  is  the  dimensionless  perturbation  to  the 
hydrostatic  pressure  p0. 

Before  deriving  the  nondimensional  equations  of  motion  and 
their  associated  boundary  conditions,  we  first  introduce  the 
following  nondimensional  variables: 
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X    = 

X 

L 

V    = 

V 

~U 

t  = 

t 
T 

H 
~H 

1 

y-l  • 

z    =    -? 

H 

-1  ■ 

w  =   _ 

'         H 

u  = 

u 

Po    = 

~Po 

P\ 

f  = 

f 

(2.19) 


where  L,  H,  U,  W,  p,,  and  T  are  the  characteristic  dimensional 
scales  defined  in  section  B  and  h  is  a  characteristic 
amplitude  of  the  topography. 

We  begin  with  the  mass  conservation  equation  (2.1)  .  Using 
(2.19)  and  remembering  that  W  =  UH/L,  (2.1)  becomes 


Ur   +  V„  +  w,   =  0  .  (2.20) 


To  assist  in  deriving  the  nondimensional  form  of  the  x- 
component  of  the  momentum  equations,  given  by  (2.2),  we  shall 
now  derive  the  nondimensional  form  for  the  convective 
derivative  D/Dt : 


JL  =  J-    +  (V+u)JL    +vi  +wJL.  (2.21 

DT         d"E         dx  by  dz 


Using  (2.19),  (2.21)  becomes 
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-5.  .1   a  +  »(V+u)  3  ♦  -?v»  ♦  £w»  (2.22) 

DT   T  Tt    L      TTx    L  ^Jy    L  TTz 


but  T  =  L/U,  therefore  (2.22)  simplifies  to 


D  ID 

~DZ         T  Dt 


(2.23) 


where 


_EL  =   ^  +  (V+U)  *  +  v  *  +  *JL.  (2.24) 

Dt    dt         dx     dy     dz 


Substituting  (2.17),  (2.18),  and  (2.23)  into  (2.2)  and  using 
(2  .19)  ,  we  obtain 


,   yS   U  Du         P\9H$ 

^>-f)-rl  +  -^T"^  =   or' 


!2.25) 


P.<«-r>£+a-o 


A  similar  analysis  for  the  y- component  of  the  momentum 
equations  (2.3)  yields 
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P0(z-{)^   +  qy  =   0.  (2.26) 


To  derive  the  nondimensional  equation  for  the  z- component  of 
the  momentum  equations  (2.4),  we  begin  with 


Pjz-n^   +  P0(z-$)9   +  ft  =  0.  (2.27) 


Substituting    (2.17),     (2.18),    and    (2.23)    into    (2.27)    and  using 
(2 .19) ,    we   obtain 


Po(z~n™^  +  P„(^-r)sr  +  gpoz  +  gri3gz  =  0.  (2.28) 


Using   T   =    L/N,H   and  W   =   N,H2/L,     (2.28)    simplifies    to 

N*ll^Jp0(z-?)V£  +  p0(z-t)9  +  9P„  +  gr/Sqr,  =  0.  (2.29) 

Using  (2.12)  and  the  definition  of  0,  0  =  N,2H/g,   (2.29) 
simplifies  to 
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£2p"(2"r)-§£  +  |(p#(«-r)-  p.(»))  +  <?z  =  0,     (2.30 


where 


e  =  — .  (2.31) 

L 


Lastly,  we  derive  the  nondimensional  vertical  velocity 
from  its  definition  (2.5).   Using  (2.19)  yields 


w  =   El.  (2.32 

Dt 


We  now  consider  the  kinematic  bottom  boundary  condition 


f  =  f         at         z  =  f.  (2.33) 


Using  (2.19),  we  obtain 


f  =  Off    at    z  =  otf,  (2.34) 


where 
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a  =  4-  (2.35) 

H 


For  the  free  surface  kinematic  boundary  condition, 

7  =  ~rj  at  z  =  H  +  rj.  (2.36) 

We  use  (2.19)  to  obtain 

£  =   -q  at  z   =  1  +  rj.  (2  .37) 

Next,  we  nondimensionalize  the  free  surface  dynamic  boundary 
condition  in  (2.9b): 

po(H+y)    +  Pp(H+y)    =po(H)       at      z  =  H  +  ^ .  (2.38) 

Using  (2.19)  and  (2.18),  (2.38)  becomes 

p0(l  +  i?)  +j8g  =  p0(l)   at   z  =  1  +  17.  (2.39) 

To  recapitulate  our  nondimensional  versions  of  equations 
(2.1-2.6)  and  the  surface  boundary  conditions  (2. 8-2. 9b),  the 
nondimensional  equations  of  motion  are 
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ux  +  vy  +  w.  =  0 ,  (2.40) 


Po(z-r)T^  +  *  =  °'  (2'41) 


Po(z~n^   +  gy  =  0,  (2.42) 


e2P.(2-n-j^  +   ^(z-n-p0(z))   +  g:  =  0,  (2.43 


W=   .££,  (2.44) 


where    £    =    H/L  and 


£t  =  -k+™4-X*v-$y+w-k- 


with  boundary  conditions 

£  =  af         at    z  =  af,  (2.46) 

f  =  17     at    z  =  1  +  r) ,  (2. Ala) 

pjl  +  rj)  +j3g  =  po(l)     at    z  =  1  +  tj  ,  (2.47Jb) 

where  a  =  h/H. 

An  exact  analytical  solution  of  these  nonlinear  partial 
differential  equations  is  not  available.   However,  we  shall 

17 


obtain  approximate  solutions  for  the  case  when  e  defined  by 
(2.31)  and  a  defined  by  (2.35)  are  both  small  positive 
constants.  Since  e  =  H/L  =  N/N,,  the  condition  e  <<  1  is 
satisfied  when  the  density  changes  rapidly  across  a  thin 
thermocline.  However,  there  is  some  ambiguity  in  the 
definition  of  the  thermocline.  We  will  consider  the 
thermocline  to  be  the  region  where  density  variations  are 
determined  chiefly  by  temperature  variations.  This  means  that 
the  submarine  is  usually  submerged  to  a  depth  below  the 
thermocline.  Is  this  a  valid  assumption?  To  avoid  acoustic 
detection  the  submarine  submerges  to  depths  below  the  mixed 
layer  or  layer  depth.  Depending  on  the  thickness  of  the 
thermocline  the  submarine  can  be  either  in  the  thermocline  or 
below  it.  We  shall  consider  only  the  case  for  which  the 
submarine  is  below  the  thermocline  in  which  the  parameter  a, 
the  ratio  of  the  characteristic  amplitude  of  the  topography  to 
the  submarine's  depth  defined  by  (2.35)  ,  is  also  small  in  this 
situation. 
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III.   THE  ASYMPTOTIC  EXPANSION  OF  THE  DEPENDENT  VARIABLES 

To  study  the  behavior  of  the  internal  waves  excited  by  a 
small  amplitude  disturbance  in  the  presence  of  a  uniform 
oncoming  flow,  as  when  a  submarine  moves  through  the  water  at 
constant  velocity  V,  we  perform  an  asymptotic  expansion  on  the 
variables  f,  q,  u,  v,  and  77  in  terms  of  the  small  parameter  a. 
This  is  reasonable  since  the  small  parameter  a  enters  the 
bottom  boundary  condition  in  an  analytic  way  and  since  the 
solutions  arise  by  analytic  processes,  it  is  natural  to  expect 
the  solution  to  exhibit  an  analytic  dependence  on  the 
parameter  ex.  Since  the  disturbance  is  0(a),   we  shall 

initially  suppose  that  the  response  is  also  0(a?)  .   Therefore, 
we  introduce  the  following  asymptotic  expansions 

f  =  a{0(x,y,z,  t)  +  a2f,  + 

q  =  cxqo(x,y,  z,t)  +  oi2qi  + 

u  =  cxu0(x,y,  z,t)  +  a2u,  +   .  .  .  (3.1) 

v  =  cxv0(x,y,  z,  t)  +  cx2v1  + 

r]  =  cx7]0{x,y,  z,  t)  +  Qf2i7,  + 

in  which  we  consider  the  limit  process  as  ex  -*   0 . 

Substituting  (3.1)  into  the  continuity  equation  (2.40)  ,  we 
obtain 
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J^iotuJ    +  -^[avj  +  ^[^(«r.)]  +  0(<*2)  +  ...  =  0.     (3.2) 


Dividing  by  a   and  taking  the  limit  as  a   -»  0,  (3.2)  simplifies 
to 


where 


and 


u  +  v  +  w  =  0 ,  (3.3) 

at      oy      oz       '  v      ' 


Q*Q  /O  A) 


_^_  =  _jL  +  V  ^  .  (3.5) 

Dct   TE    ^x 


To  show  that  D0/D0t  is  the  0(1)  part  of  the  convective 
derivative  given  in  (2.45),  we  collect  only  the  0(1)  terms  in 
the  expansion  for  D/Dt  given  in  (3.6)  below: 


+  (V  +  otu)  J-   +   av„JL  +  0(oi2)    +  .  .  .     (3.6) 


Dt       Ti  °   "53c    °~8y 
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Substituting   (3.1)   into  the  x- component  of  the  momentum 
equations  (2.41),  we  obtain 


p0(z-{)  _1  (0iuo)    +  Oiqm   +  0(cx2)    +  .  ..  =  0.         (3.7) 


Expanding  p0(z-£)  in  a  Taylor  Series  and  using  (3.1)  yields 


P„(z-f)  =  Po(z)    -  p[,{z)   [afj  +  0(a2)  + (3.8) 


Substituting  (3.8)  into  (3.7),  we  obtain 


apjz)  -J^l   +  ag^  +  0(a2)    +  .  .  .  =  0.  (3.9) 


Dividing  by  a  and  taking  the  limit  as  a  -»  0 ,  (3.9)  becomes 


Po^)^     +  Vat    =     0.  (3.10) 

0 


A  similar  analysis  for  the  y- component  of  the  momentum 
equations  (2.42)  yields 
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P°{Z)1TE  +  q<*  =   °-  (3'11} 


For  the  z-component  of  the  momentum  equations  (2.43),  we 
substitute  (3.1)  into  (2.43)  to  obtain 


t-j2  v. 

ae2^h   +  i[p„(z-r)  -  P.(*)l  +  ag^  +  0(a2)  +  ...  =  0.     (3.12) 


Expanding  p0(z-£)   in  a  Taylor  Series  as   in   (3.8)   and 
substituting  into  (3.12),  we  obtain 


D2t 
ote2-l±l        *p'0{z)S0  +  cxg^  +  0(a2)    +   ...    =  0.  (3.13) 

2/ _  U  H 


Dividing  by  a.   and  taking  the  limit  as  a  -»  0,  (3.13)  simplifies 
to 


-lp'jz)L   +  g,  =  0.  (3.14) 


Notice  that  the  £2  term  has  been  neglected.  As  previously 
discussed,  in  order  for  our  theory  to  work  we  require  £  «  1. 
Therefore,  we  see  the  £2  term  is  small.   The  fact  that  the  e2 
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term  is  small  by  itself  is  not  justification  in  dropping  this 
term  from  (3.13)  .  In  order  to  justify  the  dropping  of  the  e2 
term,  we  shall  assume  that  oie2  =  0(01)  as  a  -*  0 ,  or  that  e2  = 
o(l)  as  a  -*  0 .  Now  as  one  takes  the  limit  oe  -*  0  in  (3.13)  ,  we 
clearly  see  that  e2  -*  0 .  In  continuing  our  derivation  of  the 
z- component  of  the  momentum  equations,  we  substitute  (2.15)  in 
for  dp0/dz,  so   that  (3.14)  becomes 

g«  +  Po^Ho  =  o.  (3.15) 

Next,  we  consider  the  surface  boundary  conditions.  We 
begin  by  substituting  (3.1)  into  the  kinematic  bottom  boundary 
condition  (2.46)  to  obtain 

af„  =  oif  +  0(oi2)    +  .  .  .    at    z  =  af.  (3.16) 

Taking  the  limit  as  a  -*   0 ,  we  find 

t    =  f         at  z  =   0.  (3.17) 

Note  that  the  0(a)  bottom  boundary  condition  is  applied  at  z 
=  0.  Substituting  (3.1)  into  the  free  surface  kinematic 
boundary  condition  (2.47a)  results  in 
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ato  =  aV0  +  0{oi2)    +    .  .  .  at  z  =  1   +  0t-qo.  (3  .  18) 

Taking   the   limit   as   a.  -*  0,     (3.18)    becomes 

r0  =  Vo         at  z  =  I-  (3.19) 

Again  the  0(a)  free  surface  condition  is  applied  at  z  =  1 . 

The  last  boundary  condition  to  consider  is  the  free 
surface  dynamic  boundary  condition  (2.47b).  Substituting 
(3.1)  into  (2.47b) ,  we  obtain 

p„(l+i?)  +  a0qo   +  0(a2)    +  .  .  .  =  p0(l)   .  (3.20) 

Expanding  p0(l  +  t;)  in  a  Taylor  Series  and  using  (3.1)  gives 

pa(l+y)    =  pJD    +  ap!>(z)rJo   +  0(a2)  +...        (3.21) 
Substituting  (3.21)  into  (3.20),  we  find 

Po(l)    +  ap'0(l)ri0   +  aj8g0  +  0(a2)    +  .  ..  =  p0(l)  .      (3.22) 
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Canceling  the  p0(l)  terms,  dividing  by  a,    and  taking  the  limit 
as  a   -*  0,  (3.22)  becomes 

p,0(z)r,0   +  j8go  =  0     at     z  =  1  (3.23) 


but  using  (2.12)  simplifies  the  free  surface  dynamic  boundary 
condition  to 

0go  =  P0Vo         at  z  =  1.  (3.24) 

In  summary,  we  see  from  substituting  the  asymptotic 
expansion  (3.1)  into  (2.40-2.45),  (2.46),  and  (2.47a,b),  that 
the  0 (a)  problem  is: 


u+v+w=0,  (3.25 

ox  oy  oz  '  N 


P"(Z)TTF  +  g-  =  °'  (3-26) 


P.<*>TJlf  +  <?„  =  0,  (3.27 


<3T«  +  P^f*  =  0,  (3.28) 
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».-§£.  '  <3-29> 


where 


+  vJL,  (3.30) 


D0t       ~3~t         ~3x' 


with  boundary  conditions 

£,  =  f          at  z  =  0,  (3.31) 

ro  =  r?0          at  z  =  1,  (3.32a) 

0go  =  p0r?o          at  z  =  1.  (3.32jb) 

Combining   (3.32a)   and   (3.32b)  gives  the  single  boundary 
condition  at  the  free  surface: 

Pq0   =  pj0         at  z   =  1.  (3.33) 
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IV.   DEVELOPMENT  OF  THE  TWO-DIMENSIONAL  WAVE  EQUATION 

A.   EIGENVALUE  PROBLEM 

In  order  to  solve  (3.25-28),  (3.31),  and  (3.33),  we  look 
for  solutions  of  the  form 

u0   =  Aei(KX  +  LY'UT)  <t>'  (z) 

Vo    =  Bei(KX  +  LY-QT)    ^/(z) 

(4.1) 
r0  =  Dei{KX*LY-QT)  <f>(z) 

q0  =  Dc2eiUa  +  LY-ar)  Po  tf  (z) 

where  k  =  (K,L)  is  the  horizontal  wave  vector,  Q  is  the 
frequency,  and  c  is  the  long -wave  phase  speed  relative  to  the 
basic  state  at  rest.  Substituting  (4.1)  into  (3.25),  (3.26), 
and  (3.27)  and  solving  the  system  of  linear  equations,  we  find 
the  dispersion  relationship  to  be 

Q  =  c\k\  where  \  k  \  *  0 .         (4.2) 

Furthermore,   substituting   (4.1)   into   (3.28),   (3.31),   and 
(3.33)  results  in 
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(po4>'  (*))'  +  -^-0(z)  =   0,  (4.3a) 


with  boundary  conditions 

<t>(z)    =  0    at    z  =  0,  (4.3Jb) 

0(z)  =j8cV(z)    at    z  =  1.  (4.3c) 

The  equations  (4.3a-c)  formally  pose  an  eigenvalue  problem 
where  c'2   is  the  eigenvalue  and  </>  the  eigenf unction.  Note  that 
the  dispersion  relationship  in  (4.2)  was  obtained  for  wave 
motion  in  a  stratified  fluid  at  rest  where  V  =  0  and  D0/D0t  = 
d/dt.    In  the  problem  we  are  considering,  the  dispersion 
relationship  would  differ  from  that  found  in  (4.2)  but  the 
eigenvalue  problem  would  remain  the  same.   Hence,  our  purpose 
here  is  to  show  the  existence  and  character  of  the  solutions 
that  define  the  eigenvalue  problem.   The  solutions  have  the 
horizontal  character  of  plane  waves  travelling  in  the  (x,y) 
plane  with  a  vertical  structure  defined  by  the  eigenf unctions, 
0(z).   We  shall  see  that  the  general  solution  to  (3.25-3.28), 
(3.31),  (3.33)  can  be  obtained  in  terms  of  a  superposition  of 
normal  modes  whose  depth  dependence  is  determined  by  the 
</>s(z)  ,  the  eigenf  unctions  of  the  following  eigenvalue  problem: 
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lPo<t>«   >'  +   ^rts  =  0,  (4.4a) 

c2 


(f>s(z)    =0    at         z   =  0,  (4.4Jb) 

0,(z)  =  0c2<£K(z)    at    z  =  l.  (4.4c) 

Here  cs~2  is  the  eigenvalue,  where  cs  is  the  long- wave  phase 
speed  relative  to  the  basic  state  at  rest. 

B.   ORTHOGONALITY 

The   eigenf unctions   0s(z)   satisfy   the   orthogonality 
condition 


c^lQP04>szKdz  =  6sr  l„  (4.5) 


where  Is  is  defined  by  the  left  side  of  (4.5)  for  r  =  s.  To 
prove  (4.5)  ,  we  let  Xs  and  Xr  be  eigenvalues  with  corresponding 
eigenf  unctions  <ps  and  </>r  where  Xs  =  cs"2  and  Xr  =  cr2.  The 
differential  equations  satisfied  by  the  eigenfunctions  are 


(PA)'  +  hPo&h  =   0,  (4.6a) 

(P.**)'  +  \PoN2<t>r   =  0.  (4.6b) 
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In  addition,  both  <£s  and  $r  satisfy  the  same  set  of  homogeneous 
boundary  conditions  defined  in  (4.4b,c).  Multiplying  (4.6a) 
by  \r<pT   and  (4.6b)  by  \<ps   yields 

K(Po<t>szy<t>r    +    \\po&4>,4>r    =    0,  (4.7a) 

MP.**)'*,   +  KK  PoN24>s4>r  =0-  (4.7b) 

Subtracting    (4.7b)    from    (4.7a)    results    in 

KiPoWtr    ~    KtPo**)'*,    =     0.  (4.8) 

Through  application  of  the  product  rule 

KlPo4>*)4>r  =  M(pAf)'  -p.*«*„].  <4-9> 

(4.8)    becomes 

KiPoKW  -  \poKK  -  MpA*J'  +  \PoKK  =  o.  (4.10) 

Integrating    (4.10)    yields 


30 


(X,  -  Xr) 


r\  .  )  } 


Po4>a  <pn  dz   =  \  Po<t>n  0,    I     -  \rPo<t>sz  <pr\.  (4.11) 


Evaluating  the  right  hand  side  of  (4.11)  and  remembering  <£s(l) 
=  jffA/'^d),  0r(l)  =  /3\rl4>st(l)  ,  and  0,(0)  =  0,(0)  =  0,  (4.11) 
becomes 


(\  -  Xr)  V  p0(pa<pa  dz  =  0.  (4.12) 


If   Xs  ^  Xr,    it    follows    that 


^P04>sz<Pn  dz  =  0.  (4.13 


In  other  words,  the  z  derivatives  of  the  eigenf unctions  (0r 
and  0S)  corresponding  to  different  eigenvalues  (Xr  and  Xs)  are 
orthogonal  with  weight  p0.  Therefore,  the  normal  modes  do 
satisfy  the  orthogonality  condition  in  (4.5). 

C.   TWO-DIMENSIONAL  WAVE  EQUATION 

Since  the  eigenf  unctions,  0s(z),  form  a  complete  set,  we 
can  represent  u0,  vof  fof  and  q0  by  a  generalized  Fourier  series 
expansion.  From  (4.1)  we  know  the  general  form  of  the 
solutions  and  hence  we  put 
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Co  =  Y.As(x>y>  t)^(z)  ,  (4.14a) 


<7o  =   P0U)]Tcs2£sU,y,  t)(J)S2(z)  ,  (4.14^) 


u0  =  -J^Ds(x,y,t)4>sz(z)  ,  (4.14c) 


vQ  =  -^Es{x,y,  t)$sz(z)  .  (4.14d) 


To  derive  a  relationship  between  As  and  Bs;  we  begin  with 


C0  =  £AsU,y,  t)<t>s(z)   .  (4.15) 


We  take  the  derivative  of  (4.15)  with  respect  to  z,  multiply 
by  pocs20rz(z)  and  use  the  orthogonality  condition  in  (4.5)  to 
obtain 


I SAS  =  clC9Jloz$Bzdz.  (4.16) 

j  o 


Integrating  (4.16)  by  parts  then  gives 
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Is*s   =   C-S2P0C04>SZ|         C2J\0{9Jksz)i  dz.  (4.17) 

o  Jo 


Using    (4.4a),     (4.17)    becomes 


I        n 


ISAS  =   Cs2po(|>sz|    +    Cp0C0N2^8dz.  (4.18) 

o        jo 

Now,     we    multiply     (4.14b)     by    0re    and    from    the    orthogonality 
condition    (4.5),    we   obtain 


IS5S  =    Cdobszdz.  (4.19) 

Jo 


Integrating  (4.19)  by  parts,  we  obtain 


I       n 


i.bb  =  g04>sl   -  f  Qozbsdz-  (4.20) 

o        Jo 


Using    (3.28),     (4.20)    becomes 


I         n 


ISBS  =  qo4>s\    +    [   p0C0N2<i>sdz.  (4.21) 

o        Jo 


Combining    (4.18)    and    (4.21)    yields 
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1  1 

I 

0 


CsPotcfiszl    +   IsBs   -   flrAl    =   JsAs-  (4.22) 


Substituting  (3.33)  and  (4.4c)  for  q0(l)  and0s(l)  respectively 
and  remembering  that  fo(0)  =  f(x,y,t),  (4.22)  simplifies  to 


Bs  =  As   +  Fs,  (4.23) 


where 


ISFS  =  c*po(0)<bsz(0)f(x,y,  t)  .  (4.24) 


Note  that  Fs(x,y,t)  is  proportional  to  the  bottom  topography 
f  (x,y,t)  . 

Next,   we  substitute   (4.14a,c,d)   into  the  continuity 
equation  (3.25)  to  obtain 


-±*E^       -±***i       +  f^4>   =0.         (4.25) 


Multiplying     (4.25)    by    cs2/0o#ra    and    applying    the    orthogonality 
condition    (4.5)    then  gives 
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DQ  t  dx  dy 


(4.26) 


where  D0/D0t  is  given  in  (3.30). 

To  derive  the  two-dimensional  wave  equation,  we  use  the  x 
and  y  components  of  the  momentum  equations  (3.26)  and  (3.27) . 
Substituting  (4.14b,c)  and  (4.14b,d)  into  (3.26)  and  (3.27) 
respectively,  we  obtain 


Po 


D 


D0t 


E  w* 


+    Po-T- 
OX 


E  c'  bs  K 


=  o, 


(4.27a) 


P<> 


Dt 


E  Es  ** 


+  p 


-by 


HC'Bs*s 


{A. 21b) 


Multiplying  (4.27a)  and  (4.27b)  by  cs20„,  and  applying  the 
orthogonality  condition  (4.5)  then  yields 


D0DS  2  dBs 


0, 


4.28a) 


D0ES  2  dBs 

-DT        C'Sy 


=  0 


:4.28Jb) 
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We  next  take  the  x  derivative  of  (4.28a)  and  the  y  derivative 
of  (4.28b)  and  add  the  two  results  to  obtain 


Dt 


dDs         dEs 

si  + 


-57 


+  c2rV2B„   =  0 


(4.29) 


Substituting  (4.23)  and  (4.26)  into  (4.29)  finally  yields 


fa, 


D, 


p  ■  *>'[*.  +  '.]• 


(4.30) 


where 


and 


V2  = 


dx7 


dy2 


:4.3i: 


Di 


Dt' 


4-  +v* 
Tt    ^x 


12 


(4.32) 


Equation  (4.30)  is  the  inhomogeneous  two-dimensional  wave 
equation  whose  solution  will  be  obtained  in  Chapter  V. 

To  summarize,  we  see  that  substituting  (4.14a-d)  into 
(3.25-3.27)  implies 


36 


B,  =  A,  +  F„  (4.33a! 


AA         3D5    +    dE, 
D0t        ~3x        Ify  ' 


4.33Jb) 


_f_i    =   cfV2[  A,   +  F,  ]  ,  (4.33c) 

D  X2 


where 


ISFS  =  c2Po(0)<t>sz(0)f(x,y,t)  ,  (4.33d) 


J?L    =  ^L    +  V*  (4.33e) 

D,t         dt  dx 


and 


V2  =   -2-    +   _^_.  (4.33f) 

3x2         dy2 


Here,    since   f   vanishes   at    t    =    0,    Fs(x,y,0)    =    0   and   from    (2.7) 
we   have 

Fs(x,y,t)    ~  Gs(x,y)  as  t  -»  oo,  (4.34a) 

where 
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ISGS   =  c^po(0)<t>sz(0)g(x,y)   .        (4.34jb) 


The  initial  conditions  for  (4.14a-d)  are  that  A,.,  Bs/  Ds,  and 
E„  all  vanish  at  t  =  0 . 
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V.   SOLUTION  TO  THE  TWO-DIMENSIONAL  WAVE  EQUATION 

A.   METHOD  OF  DESCENT 

To  solve  (4.33c)  we  first  introduce  the  following  change 
of  variables: 


f  =  x  -  Vt, 

y  =  y,  (5.ia) 

a   =  t. 


Derivatives  with  respect  to  x,  y,  t  can  be  expressed  as 
derivatives  in  the  new  variables  f,  y,  o   as  follows: 


a       a 
^x  sr 

d     =      d 
~&y       ~&y' 

a  d  d 

Tt   To  Tf 


(5.1b) 


Substituting  (5.1a,b)  into  (4.33c)  yields 


a2A        2    a2A,      a2Av , 
a  o-2  a  ^2      ay2 
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where 


f.it+Vo.y.a)    =  c]  [  JL   +  JL  ]  F.if+Vo.y,  a)  . 


(5.3) 


We  shall  solve  (5.2)  by  beginning  with  the  solution  to  the 
following  three-dimensional  wave  equation: 


d2u        _2  r    d2u         d2u         d2u 


dt' 


=  c2[ 


dx2        dy2        dz: 


], 


u(x,y,z,0)    =*(x,y,z), 


ut(x,y,  z,0)    ='*,(x,y;z). 


(5.4a) 

(5.4Jb) 

(5.4c) 


Using   Kirchoff's    formula    (see   Strauss    1992),    the   solution   to 
(5.4a-c)    is 


u(x0,y0,z0,t0)    =  L—   [7*(x,y,z)dS 

47TC2t.   j,J 


a 


:7TC2t0  4,J 


(5.5) 


,y,z)dS 
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where 


2  «-2 


D  =  [  (x,y,z)   :  (x-xj2  +  (y~yo)2   +  (z-z„)2  =  c 2t„  ] 


(5.6) 


Note  that  the  value  of  u  (x0,y0,  z0,  t0)  depends  only  on  the  values 
of  SMx,y,z)  and  <t>(x,y,z)  on  the  spherical  surface  D. 

By  the  method  of  descent,  we  now  obtain  the  solution  to 
the  following  two-dimensional  wave  equation: 


d2u        „2  r  d2u         d2u 

=  C     [   +  

dt2  dx2         by2 


]  , 


u{x,y,  0)  =  $(x,y)  , 


(5.7a) 


(5.7Jb) 


ut(x,y,0)    =  *(x,y)  . 


5.7cl 


We  regard  the  solution  u(x,y,t)  of  (5.7a-c)  as  a  solution  to 
the  three-dimensional  problem  (5.4a-c)  with  initial  data 
independent  of  z.   By  Kirchoff's  formula  (5.7),  we  have 


u(0,0,  t) 


47TC2t  1 


*(x,y)dS 


(5.8) 


-Jt 


_L_  f  \${x,y)dS 

:7TC2t  JJ 
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where 


Dt  =    [  (x,y,z)   :  x2  +  y2  +  z2  =  c2t2]   .      (5.9) 

This  is  the  solution  to  (5.7a-c)  but  it  can  be  simplified.  To 
simplify  (5.8),  we  must  express  the  surface  integral  as  an 
ordinary  double  integral.  It  is  twice  the  integral  over  the 
top  hemisphere  z  =  (c2t2-x2-y2) m  where  (x,y)  ranges  over  the 
region  Dt  =  {  (x,y)  :  x2  +  y2  <,   c2t2  }.   By  definition 

dS  =    | Rx  x  R   | dxdy ,  (5.10a) 


where 


\x#y/V/(c2t2-x2-y2)  /,  (5.10Jb) 


and 


Rx  x  Ry   |  =  ^ (5.10c) 

\/(c2t2-x2-y2) 


Substituting  (5.10a-c)  into  (5.8)  yields 
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u(0,0; t)    = 


27TC 


*{x>y]    dbcdv 


Tt 


d?  yc2t2-x2-y2 


$(x,y) 


2ttc 


'f 

d f  \lc2t2-x2-y2 


dxdy 


(5.11) 


This  is  the  solution  at  the  point  (0,0,t).   Therefore,  the 
solution  at  a  general  point  is 


u(x0,y0,t0)    -   JLJ 


*(x,y) 


kJ  jc2tl-{x-x0)2-{y-yo)2 


dxdy 


TE~ 


—  ff 

!7TC  J  J 


<J>(x,y) 


*  J  /c^-U-x^My-y,)2 


dxdy 


5.12 


where 


D,  =  [  (x,y)  :  (x-x„)2  +  (y-yj2   s  c2t2].       (5.13) 

Note  that  (5.12)  shows  that  the  value  u(x0,y0/t0)  depends  on 
the  values  of  $(x,y)  and  ^(x,y)  inside  the  disc  defined  in 
(5.13).  This  illustrates  the  well  known  fact  that  Huygens ■ 
principle  fails  in  two  dimensions.  We  shall  use  (5.12)  in  the 
solution  to  the  nonhomogeneous  two-dimensional  wave  equation 
discussed  in  the  next  section. 
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B.   DUHAMEL'S  PRINCIPLE 

We   now   consider   the   following   nonhomogeneous   two 
dimensional  wave  equation 


VU  c2[dhl       +        dhj  =g{Xfyit)i  (5>14a) 


dt2  dx2        dy2 

with  homogeneous  initial  conditions 

u(x,y,0)  =  0,  (5.14Jb) 

ut(x,y,0)    =   0.  (5.14c) 

By  Duhamel '  s  principle,  the  solution  to  (5.14a-c)  is  given  by 


u(x,y,  t) 


U(x,y,  t,  s)  ds  (5.15 


where    for   each   s    &    0,    U(x,y,t,s)    is    the   solution   of 


d2u      c2[£u  +    Pu}   __0       for    ^  (5>16a) 


dt2  dx2         dy2 


with  initial  data  prescribed  on  the  plane  t  =  s  such  that 
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U{x,y,  s,  s)    =   0,  (5.16Jb) 

U,[x,y,s,s)    =g{x,y,s).  (5.16c) 

Since  the  wave  operator  is  invariant  under  translations  in 
time  t  (see  Strauss  1992),  we  can  shift  the  time  t  by  a 
constant  and  still  have  a  solution.   We  let  V(x,y,t,s) 
U (x, y, t+s , s) ,  so  that  V(x,y,t,s)  is  a  solution  to 


Va    ~    c2<V*r  +  Vyy)      =    °     f°T  t  i  0,  (5.17a) 


with  initial  data  prescribed  on  the  plane  t  =  0  such  that 

V(x,y,0,s)    =U(x,y,s,s)    =0,  (5.17Jb) 

Vt(x,y,0,s)    =  Ut{x,y,s,s)    =  gr(x,y,  s).  (5.17c) 

Using     Kirchoff's     formula      (5.12)      to     obtain     a  solution     to 
(5.l7a-c)    then  yields 

V{x,y,t,s)    =  -iL  f  f              gr(x,y,  s)             _  & d~  (5.18) 

27TC  J    J      l„7t.7_  #  —  _vn  2_  /  — _„\  2 


b/  \jc2t2-  (x-x)2-  (y-y)2 
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where  D,  =  {  (x,y)  :  (x-x)2  +  (y-y)2  <;  c2t2  }.  Since  our  purpose 
here  is  to  solve  the  nonhomogeneous  two-dimensional  wave 
equation  in  (5.14a-c)  by  Duhamel ' s  principle  (5.15),  we  can 
again  shift  the  time  t  by  a  constant  so  that  U(x,y,t,s)  = 
V(x,y, t-s, s) .  This  last  expression  for  U(x,y,t,s)  together 
with  (5.18)  then  implies 


U(x,y,  t,  s) 


lire  J  J 


g(x,y,s) 


]d  }  \jc2  ( t-s)2-  (x-x) 2-  (y-y) 2 


dxdy 


5.19 


where 


D,  =    [    (x,y)    :  (x-x)2  +  (y-y)2   s  c2(t-s)2  ] 


(5.20) 


Now  using  Duhamel ' s  principle  (5.15),  we  obtain 


u(x,y,  t)    = 


2ttc 


g(x,y,s) 


d)   V 'c2  (t-s)2-  (x-x)2-  (y-y)2 


dxdy 


ds     (5. 21' 


the  solution  to  (5.14a-c)  where  D,  is  defined  in  (5.20). 

In  deriving  the  solution  to  the  nonhomogeneous  two- 
dimensional  wave  equation (5 . 14a-c) ,  we  have  solved  (5.2).   We 
change  the  names  of  the  dependent  and  independent  variables  in 
(5.21)  to  obtain  the  solution  of  (5.2): 
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A,(f,y,ff)  = 


27TC. 


fs(x+Vs,y,  8) 


°.    \lcl(a-s)2-(x-{)2-(y-y)2 


dxdy 


ds         (5.22) 


where  D„  =  {  (x,y)  :  (x-f)2  +  (y-y)2  s  c2(a-s)2  }.   Letting  f  = 
x  -  Vt  where  o   =  t,  (5.22)  becomes 


As(x,y,  t)    = 


27TC. 


[{ 


fs(x+Vs,y,  s) 


fiTi 


dxdy 


t-s)2-(x-(x-Vt)  )2-{y-y) 


ds 


5.23 


where  Dt  =  {  (x,y)  :  (x-(x-Vt))2  +  (y-y)2  s  cs2(t-s)2  }. 

To  derive  the  solution  for  the  nonhomogeneous  two- 
dimensional  wave  equation,  we  began  by  solving  the  homogeneous 
three-dimensional  wave  equation  using  Kirchoff's  formula 
(5.5)  .  Next,  we  used  the  method  of  descent  to  solve  the 
homogeneous  two-dimensional  wave  equation  by  reducing 
Kirchoff's  surface  integral  formula  to  an  ordinary  double 
integral  (5.12)  .  We  developed  the  solution  to  the  homogeneous 
two-dimensional  wave  equation  in  order  to  solve  the 
nonhomogeneous  two-dimensional  wave  equation  by  applying 
Duhamel ■ s  principle  (5.21).  This  gave  the  solution  of  the 
wave  equation  in  (5.2)  in  the  form  given  above  in  (5.23) .   In 
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Chapter  VI  we  will  examine  the  solution  to  (5.23)  for  long 
time. 
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VI.  LONG  TIME  SOLUTION 

A.   Cs2  *   V2 

Before  looking  at  the  behavior  of  (5.23)  as  t  -»  oo  (where 
the  topography  fs(x+Vs,y,s)  has  a  general  form),  we  shall 
first  introduce  a  specific  moving  point  source  for  the 
topography,  so  that 

fs(x,y,s)    =  6(x-x0)  6(y-y0)  g(s)  .  (6.1a) 

Consequently 

fs(x+Vs,y,s)    =  6(x+Vs-x0)  6(y-y0)  g(s)  ,        (6.lJb) 

where  g  is  a  non-negative  function  satisfying 

gr(0)  =0    and         g(s)    =1    for         s   a  s0.  (6.1c) 

If  we  can  solve  (5.23)  by  substituting  (6.1b)  in  for 
fs  (x+Vs,  y,  s)  ,  then  we  can  solve  (5.23)  for  any  form  of 
topography  by  applying  the  principle  of  superposition. 

To  solve  (5.23),  we  begin  by  substituting  (6.1b)  into 
(5.23)  to  obtain 
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As{x,y,  t) 


27TC 


V/ 


6  (x+Vs-x0)  6  (y-y0)  g(s) 
c)  ( t-s) 2-  [x-  (x-Vt)  ] 2-  (y-y) 


dx  dy 


(6.2) 


ds 


We  let  X  =  (x+Vs-xJ  and  Y  =  (y-y0)  and  substitute  into  (6.2) 
to  obtain 

As{x,y,  t)  = 


27TC. 


6(X)  6(7)  gr(s)  dZdY 


JJtt— 

V  ^/c*  ( t-s) 2-  [*+V(  t-s)  -  (x-x0)  ]2-[Y-  (y-y0)  ] 2 


(6.3) 


ds 


where 


D,   = 


2    _    2 


(X,  y)  :  [x+v(t-s)-(x-x0)]2   +  [r+y0-y]2  *  c;(t-s) 


6.4 


In  order  to  solve  (6.3)  we  must  determine  those  s  for  which 
(6.4)  is  satisfied.   Letting  t  =  (t-s)  for  (t-s)  ;>  0,  (6.4) 
simplifies  to 


.2_T/2\  T2 


(c;-V2)t2  +  2V(x-x0)t   -  [(x-xj2  +  (y-y0)2]    *  0 


(6.5) 
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Note  that  X  and  Y  are  not  included  in  (6.5)  since  the 
integration  in  (6.3)  with  respect  to  X  and  Y  makes  both  X  and 
Y  become  zero.   Solving  (6.5)  for  t  yields 

-V(x-x0)   ±  sjc2s{x-x0)2   +  (c^-V2)  (y-y0)2  ^ 

T.    =  ! .         (6.6) 

C-V2 


The  values  of  t+  and  t,   will  depend  on  whether  cs2  >  V2  or  cs2  < 
V2.   We  shall  examine  both  cases  in  the  following  sections. 
1.   Cs2  >  V2 

When  the  long-wave  phase  speed  is  greater  than  the 
velocity  of  the  submarine,  we  see  from  (6.6)  that  t+  >  0  and 
t_   <  0 .   Factoring  (6.5)  yields 


,2_tt2 


(c,-V2)   (T-T  +  )  (T-T.)  ;>  0  (6.7) 


where  r+  and  r.  are  defined  in  (6.6)  .  Since  the  left  hand  side 
of  (6.7)  must  be  greater  than  or  equal  to  zero,  t  must  be 
greater  than  or  equal  to  t+  (t  2  r+)  .  Hence,  s  must  be  less 
than  or  equal  to  (t-T+).  Therefore,  the  upper  limit  of  s  in 
(6.3)  becomes  (t-r+).  Now,  integrating  (6.3)  with  respect  to 
X  and  Y  yields 
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AM.y.t)    -   X  -[•'-•- g(s)ds  (6.8) 


where  r+  and  t.  are  defined  in  (6.6).   Using  (6.1c),  we  break 
(6.8)  into  the  following  parts 

As(x,y,  t)    =   J,  +  J2,  (6.9a) 


where 


I,  = 


27TC5(C'-V2)1/2 


gr(s)  ds 


0  \l((t-T.)s)((t-T_)  ~S) 


(6.9b) 


and 


J2  - 


1 p.         ds  (6  >9c) 

2irc,<c,a-V)'"-k   y((t_T+)_s)((t_Tj_s) 


We  first  look  at  the  behavior  of  I,  as  t  -»  <».  We  see 
from  (6.9b)  that  for  long  time,  I,  is  0(l/t).  Therefore,  we 
can  neglect  this  portion  of  As(x,y,t).  Now  in  evaluating  I2 
for  long  time,  we  let  r  =  (t-s)  so  that  I2  simplifies  to 
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2*c,{cl-W*h        Mr-T.)(T-T_) 


We  complete  the  square  of  the  denominator  of  (6.10)  to  obtain 


dT  (6.11) 


2     27TCJ(C,2-l'2),/2>-    /77 


-a)z-b 


2_K2 


Here  a  =  (t++t.)/2,  and  b2  =  (t  +  -t.)2/4.   Letting  bu  =  (r-a)  and 
substituting  into  (6.11)  then  yields 


I2-    1 f     dU    .  (6.12) 


Integrating  I2,  we  find 
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^2     = 


2ttcs(c;-V) 


2\  1/2 


Ln| 


t-s-a] 


2ttcAc;-Vz) 


•2\  1/2 


in  | 


r-a 


t-s-a] 


-  1 


(T+-a) 


1/2 


1/2 


(6.13) 


For  long  time  t,  I2  is  0(ln  t) . 

In  the  case  where  cs2  >  V2,  we  see  that  the  long  time 
solution  for  A,.(x,y,t)  is  0(ln  t)  .   Next  we  examine  the  case 
for  which  cs2  <  V2. 
2.   Cs2  <  V2 

When  the  long -wave  phase  speed  is  less  than  the 
velocity  of  the  submarine,  we  see  from  (6.5)  that  x  must  be 
greater  that  x0  since  t  is  greater  than  or  equal  to  zero.  We 
must  reexamine  (6.6)  to  see  if  a  positive  real  value  for  t  is 
possible.   Rearranging  (6.6)  for  cs2  <  V2  yields 


r±  = 


V(x-x0)    ±  Jc2Ax-x0)2  -   (V2-c2s)  (y-yy 
V2-c2s 


(6.14) 
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From  (6.14)  we  see  that  a  real  solution  for  t   exists  only  if 
the  following  inequality  holds 


cUx-x0)2   -  {V2-c2s)  (y-yj2  *   0 


(6.15) 


Hence,  a  solution  exists  only  within  the  region  defined  by 


x  &  x„  + 


V2 


-  1 


1/2 


(y  -  Y0)  for         y  >  yo   , 


6.16a! 


and 


x  *  xa   + 


_2 


-  1 


1/2 


(y0  -  y)  for         y  <  ya 


(6.16Jb) 


Since  t+  and  r.  are  both  positive,  we  see  from  (6.7)  that  t 
must  be  between  t+  and  r.,  ie.,  (t.stst  +  )  .  Therefore,  s  must 
be  less  than  or  equal  to  (t-T.)  . 

To  solve  ^(x^jt)  for  cs2  <  V2  as  t  -*  oo,  we  start  with 
(6.3) .   Integrating  (6.3)  with  respect  to  X  and  Y  yields 
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AM.y.t)    -   !L_  [- g(s)dS  (6.17) 

2TTCS(V2-CS)  j>     y((t-Tj-s)((t-T.)-s) 

where  t+  and  t_   are  defined  in  (6.14)  .   Using  (6.1c),  we  break 
(6.17)  into  the  following  parts 

As(x,y,  t)  =  J3  +  J4  ,  (6.18a) 

where 


I,  = 


3    o^~  nr2_„2\  1/2 


27r^(v2-^)1/2V((t-rJ-s)((t-rJ-s) 


gr(s)ds  (6.18b) 


and 


x4  =  * '  f  ■      ds      

2ttc,(c^)"2  J  ^/((t_T  +  )-s)((t_T_)_s) 


(6.18c) 


Looking  at  the  behavior  of  I3  as  t  -*  oo,  we  see  from 
(6.18b)  that  I3  is  0(l/t).   Therefore,  we  can  neglect  this 
portion  of  As(x,y,t)  .   The  evaluation  of  I4  is  similar  to  the 
analysis  as  developed  for  I2: 
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1/2 


^4     = 


2ttcs(v'-c;) 


2\  1/2 


Ln 


t-s-a! 


(t-s-a) 


-  1 


(6.19) 


1/2 


-1 


2-ncs(Vz-cs 


2\  1/2 


In 


(r.-a) 


(r.-a) 


Here  a  =  (t  +  +  t.)/2  and  b2  =  (t+-t.)2/4.   For  long  time,  I4  is 
0(ln  t) . 

In  the  case  where  V2  >  cs2,  we  see  that  the  long  time 
solution  for  A^Xjy^t)  is  0(ln  t)  .  Next,  we  examine  the 
resonant  case  where  cs2  =  V2. 

B.   Cs2  =  V2 

Previously,  we  developed  the  long  time  solution  for  the 
cases  when  cs2  >  V2  and  cs2  <  V2  with  the  interesting  result  that 
in  both  cases  the  solution  had  growth  of  0(ln  t) .  The  reason 
for  this  growth  is  the  failure  of  Huygens ■  principle  in  two 
dimensions.  Though  ln(t)  tends  to  infinity  for  long  time,  it 
is  a  function  of  extremely  slow  growth,  and  for  the  times  we 
are  interested  in,  it  can  be  considered  as  a  bounded  function. 

In  the  resonant  case,  we  are  concerned  with  the  behavior 
of  A,.(x,y,t)  as  t  -»  oo  for  the  case  when  the  speed  of  the 
submarine  is  equal  to  the  long- wave  phase  speed  (cs2  =  V2)  . 
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To  determine  if  a  solution  to  (6.3)  is  possible  for  cs2  = 
V2,  we  must  reexamine  (6.5).  From  (6.5)  we  see  that  a 
solution  exists  only  if  s  satisfies  the  following  inequality 


(x-x  ) 2  +  (v~v  ) 2 

s   s  t  -  ^— -^-  for  V(x-x0)    >  0.      (6.20) 

2V{x-xJ 


Integrating  (6.3)  with  respect  to  X  and  Y  yields 

As{x,y,t)    = 

1  [*  g(s)  ds  (6.21 


—  f 

!7TC  Jo 


^Vix-x^)  (t-s)    -   [(x-xj2  +   (y-y0)2] 


Here  A  =  [t  -  (  (x-x0)  2+  (y-y0) 2)  /  (2V  (x-xj  )  ]  .   Using  (6.1c),  we 
break  (6.21)  into  the  following  parts 

As(x,y,  t)    =  J5  +  J6  ,  (6.22a) 


where 
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J5  = 


2  7TC.. 


g(s)  ds 


(6.22b) 


^2V(x-x0)  (t-s)        [(x-xj2  +   (y-y0)2] 


and 


J6  = 


^  J 


ds 


(6.22c) 


\J2Vix-Xo)  (t-s)    -  i(x-x0)2   +  (y-y„)2] 


Here  A  =  [t  -  (  (x-x0)  2+ (y-y0) 2)  /  (2V  (x-xj  )  ]  . 

Looking  at  the  behavior  of  I5  as  t  -»  oo,  we  see  from 
(6.22b)  that  for  long  time,  I5  is  0(l/t1/2).   Therefore  we  can 
neglect  this  portion  of  A^XjVjt)  .  Integrating  I6  directly,  we 
obtain 


It   = 


i  \J2V(x-x0)  (t-s0)   -   [(x-x0)2+  (y-y0)2] 


(6.23) 


2ttc.  V(x-x) 


From  (6.23)  we  see  that  I6  is  0(t1/2)  as  t  -*  oo.  Therefore,  we 
see  that  the  solution  for  As(x,y,t)  at  resonance  grows  like 
0(t1/2)  for  long  time. 
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To  derive  the  long  time  solution  (t  -*  oo)  for  A^x^t)  ,  we 
began  by  introducing  a  moving  point  source  for  the  topography, 
fs  (x+Vs,  y ,  s)  ,  with  the  knowledge  that  if  a  solution  was 
attainable  by  using  (6.1c)  ,  we  could  solve  (5.23)  for  any  form 
of  the  topography  by  the  principle  of  superposition.  We  first 
looked  at  the  case  where  cs2  >  V2  and  determined  that  the  long 
time  solution  for  As(x,y,t)  was  0(ln  t)  .  Next,  we  looked  at 
the  case  where  cs2  <  V2  and  saw  that  a  solution  existed  only 
within  a  certain  region  defined  by  (6.16a,b).  The  long  time 
solution  for  As(x,y,t)  in  this  case  was  also  0(ln  t)  .  We 
attribute  the  0(ln  t)  growth  to  the  failure  of  Huygens ' 
principle  in  two  dimensions.  For  the  resonant  case  where  cs2 
=  V2,  we  found  the  solution  for  Aj(x,y,t)  to  be  0(t1/2)  for  long 
time.  From  this  analysis,  we  see  in  this  problem  that  if  a 
submarine  patrols  at  or  near  a  long-wave  phase  speed  a 
significant  resonant  could  occur  for  long  times. 
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VII.   CONCLUSIONS 

This  thesis  has  presented  a  theoretical  study  of  the  flow 
of  a  stratified  fluid  over  variable  topography  as  a  model  of 
the  resonant  interaction  of  a  submarine's  wake  with  a 
stratified  fluid.  We  considered  the  case  in  which  a  submarine 
patrols  at  a  speed  near  the  long-wave  phase  speed  of  an 
internal  wave.  Internal  waves  are  generated  by  the  flow 
around  the  submarine ' s  wake . 

The  long  time  solution  obtained  in  our  model  for  the  cases 
when  cs2  >  V2  and  cs2  <  V2  had  growth  of  0(ln  t)  .  (see  (6.13)  and 
(6.19)).  We  attribute  the  0(ln  t)  growth  to  the  failure  of 
Huygens '  principle  in  two  dimensions.  For  the  resonant  case 
where  cs2  =  V2,  we  found  the  solution  to  be  0(t1/2)  for  long 
time.  (see  (6.23))  .  It  is  important  to  note  that  for  the 
time  scales  we  are  concerned  with,  both  long  time  solutions 
for  the  nonresonant  cases  are  bounded  and  have  extremely  slow 
growth.  Therefore  it  is  the  resonant  case  which  is  important. 

The  growth  in  time  occurring  in  the  resonant  case 
indicates  that  perhaps  nonlinear  effects,  balanced  by 
dispersion,  could  cause  significant  upstream  disturbances. 
From  these  disturbances,  the  possibility  exists  of  detecting 
a  submarine  by  the  internal  waves  its  wake  generates. 
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Further  studies  in  this  problem  should  begin  by  rescaling 
the  asymptotic  expansions  to  balance  the  leading  order 
quadratic  nonlinear  terms  with  the  topographic  forcing.  This 
will  achieve  a  balance  between  nonlinearity  and  dispersion. 
We  expect  the  solution  to  this  problem  to  show  that 
significant  upstream  disturbances  could  occur  when  the 
submarine  patrols  at  a  speed  at  or  near  a  long -wave  phase 
speed. 

If  internal  waves  were  generated  only  by  underwater  moving 
objects  then  the  problem  of  detecting  the  submarine  by  the 
internal  waves  its  wake  generates  would  be  made  simpler,  but 
this  is  not  the  case.  Internal  waves  are  a  natural  occurring 
phenomenon  in  the  oceans.  This  complicates  the  problem 
because  in  order  to  detect  the  submarine  by  the  internal  waves 
that  it  generates,  the  ambient  internal  wave  patterns  at  the 
location  of  the  submarine  must  be  known.  The  problem  of 
discerning  ambient  internal  waves  from  internal  waves 
generated  by  the  submarine's  wake  is  similar  to  the  problem  of 
detecting  a  submarine  passively  (such  as  with  a  tail).  The 
sonar  operator  filters  out  the  ambient  noises  or  background 
noises  of  the  ocean  from  the  incoming  signal  to  detect  the 
acoustic  signature  of  the  submarine.  From  further 
understanding  of  the  ocean's  naturally  occurring  internal 
waves,  we  may  one  day  be  able  to  discern  internal  waves 
generated  by  a  submarine  from  those  of  the  ocean. 
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